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ABSTRACT 



Many methods have been proposed and tested in the 
past to obtain an optimal response for step inputs to 
automatic feedback control systems. Most of these 
methods have led to the use of sophisticated control 
devices ranging from small analog to large digital 
computers. Here the possibility of a simplified 
switching logic combined with an open-closed loop 
servomechanism is investigated. Deadbeat response to 
step inputs was the object of this study rather than 
a time optimal response. Two types of logic were 
investigated. A time invariant controller was analyzed, 
built and tested. The system works on the principle 
of constant switching times with the output being 
controlled by an open locp driving voltage which is 
proportional to the input step size. At the completion 
of the open loop mode of operation, the system is 
returned to the normal closed loop mode. 

The writers wish to express their appreciation 
for the assistance and encouragement given them by 
Dr. John Ward, of the U. S, Naval Postgraduate School, 
in this investigation. 
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1 . Introduction, 

Numerous schemes have been applied to positioning 
feedback control mechanisms to obtain optimum response 
for step inputs. Some of the methods used have been 
(1) dual mode operation 1 using a relay control computer 
which approximates the optimal switching lines, (2) use 
of compensation networks to improve relay performance , 
(3) application of discontinuous damping to a relay 
servo^, (4) conditional switching techniques'" 1 , and 
others 5,6 c j n general* the basic idea behind these 
various approaches is the desire to obtain deadbeat 
response in minimum time for a step input. 

The time optimization requirement inherently 
requires a sophisticated controller that can compute 
the proper switching times accurately. A question 
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^Harris, McDonald, Thaler, Quasi-optimization of 
relay servos by use of discontinuous damping, Applica- 
tions and Industry, November, 1957. 
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system for space vehicles, Aerospace Engineering, Oct. 1962 

^T, R. Frederickson, A time-optimal positioning servo. 
Control Engineering, February, 1963. 

^G. J. Thaler, M. P. Pastel, Analysis and design of 
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which might he asked is: "If the requirements for 

time optimization were dropped, retaining the provision 
for deadbeat response, would it be possible to design a 
simple controller for a second order system that would 
be of practical value?" 

In this paper two approaches to the development of 
a simple controller for a second order system are 
investigated. In both cases the system operates in 
two modes, linear and nonlinear, which correspond to 
closed and open loop modes respectively. In the non- 
linear or open loop mode a driving voltage is applied to 
the plant for a prescribed amount of time. The voltage 
is then reversed until the system reaches an output 
which either equals or approximates the input. At this 
point the controller returns the system to the linear, 
closed loop, operation. 

As previously stated, in both approaches investi- 
gated, relay switching is caused to occur on a time 
basis. One of these approaches sets the switching 
times proportional to the command signal while the other 
holds the switching times constant while setting the 
driving voltage. The latter scheme is suggested by the 
principle of superposition. 

The system of Pig. 1 is the basic block diagram for 
both approaches. The step input commands are represented 
by R, the system output position by C and the output, of 
the open loop controller by V. The plant is second order 

apPao 



and has the transfer function: 



o) 



K/f 

5 ( s 'Xr) 



The system is to be controlled by an open loop controller 
which functions on command signals only. Any load 
perturbations or other similar disturbances within the 
loop will not affect the controller. 




Fig, 1 - Basic System 

The assumptions made for the theoretical investi- 
gations are (1) the system is linear, (2) the system is 
second order, (3) command signals are to be steps only 
and (4) new commands are not to be given until the 
previous command has been completely executed. 



2. Time proportional to command signal. 

With the output voltage of an open loop controller 
a constant in magnitude, there is a specific relation- 
ship for the time to reverse the polarity of the driving 
voltage and the time to return the system to normal 
linear operation in order to force deadbeat response. 

By using an approximation to this relationship, a theore- 
tical controller was derived. The system Incorporating 
this controller was then investigated on the CDC digital 
computer. 

For purposes of this investigation, the proposed 
system is to function as follows? At some time, t Q , a 
step command signal, + R, is applied to the system. 
Application of the command signal causes the loop to 
be opened and a voltage, + V, to be applied to the plant 
for a period of time, t s . At time t g the polarity of 
the signal is reversed for a period of time, t . The 
total time of nonlinear operation, t t , is defined by? 

(2) t t = t g - t r 

At time t^ the controller closes the loop returning 
the plant to standard closed loop operation. If at 
time t^, the plant output has arrived at the desired 
position commanded by the input, then true deadbeat 
response has been achieved. 

A digital program was written to compute the times 
t s and t^ for various values of inputs. The system 
parameters, K and , as well as the applied voltage, +V, 
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were held constant,. For the system under evaluation 
the following values were arbitrarily selected? 

K/f = 1, l/<f = 5 and V = 100. Each combination 
of t g and t^ that results in deadbeat response specifies 
a command signal R. Thus, from the program data, a 
plot of t 0 and t^ versus R was made for 05R^TT radians 
as shown in Fig. 2. From these data, t and t can be 
approximated by a straight line as shown in the figure. 
These straight line approximations of t g and t^ can be 
evaluated in point slope form ass 

(3) i s = .075j*\ + 0.05 

(4) t t ,092j>?l + O.IOS 

Since t g and t^ are independent of the sign of the 
input, the magnitude of R is specified in the above 
equations . 

It is immediately obvious from Fig c 2 that the 
straight line approximation is extremely poor for small 
inputs, but this problem can be eliminated if the 
system is held in closed loop for small inputs. Thus, 
in this system, the minimum step size for open loop 
operation would be about .25 radians. 

Using the switching times determined from equations 
(3) and (4), a program was written and appended to "Program 
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Fig. 2. Command signal varsue switching times 



I , 




I 

I 






i 



I 




Pig. 10. Transient response of 
time proportional to command 
controller. Input equals two 
radians • 



TIME 




I 






7 

Analog" to study the performance of a system utilizing 
this type controller. Pour runs were made with inputs 
of 1.5» and 2 radians. The resulting phase 

planes and time response characteristics as plotted by 
the computer are shown in Figures 3 to 10 . By inspection 
of Pig. 2 it is possible to predetermine whether 
switching will occur early or late. For examples if 
the input is ,5 radians 9 Pig, 2 predicts that both 
switching times* t s and t r * will be late and the system 
will overshoot. This is verified in Figures 3 and 4. 

If the input is 1.5 radians* it would be predicted 
that switching would be early. This is verified in 
Figures 7 and 8. 

This approach does offer some improvement in 
time“tO"Steady-state over normal linear operation but 
it is limited to a comparatively small range of inputs, 

A controller that could vary the times of switching as 
specified by equations (3) and (4) would not be simple 
to construct* which is contrary to one of the main 
design objectives as set forth in section 1. 

##*##*#* 

^"Program Analog" was written by Dr, J» R« Ward 
of the U, S. Naval Postgraduate School, The program as 
used is shown in Appendix I. 
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Time invariant controller 



Assume that a given linear second order system is 
controlled in such a manner as to provide a deadbeat 
response for a step Input. If the output of the linear 
system is proportional to a step input* then by the 

O 

principle of superposition * the output is doubled 
when the input is doubled. Furthermore* if the 
acceleration time* t c , * and the deceleration time t , 
are known for a specified step command* then it 
will be possible to obtain deadbeat response for all 
step commands if the switching times are held constant 
while the driving voltage is varied in proportion to 
the step command. 

To determine the appropriate switching times* it 
is first noted that if the controller voltage is V* 
then the system output will* If the Initial conditions 
are zero, be» 



’ 3 ) 



C(s) - 



kvA_ 

S 2 (S + /*) 



The step voltage +V is to be applied for a finite time, 
t g , and the step voltage “V Is applied for the finite 



time, t 



If at t s + the output of the system is 



r* — ’'s t 
identically equal to the Input* and the output rate is 

zero, then system response is Indeed deadbeat. It is 



» # « * # # 

O 

M. E. Van Valkenburg, Network Analysis* Prentlce~ 
Hall* Inc., pp. 79°80. 
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now necessary to find the required conditions for such 
deadbeat response.. 

The inverse transform of equation (5) is; 

( 6 ) CM = " 0 

To obtain the system rate at time t it is only necessary 

to take the derivative of equation (6). Thus* 

(7) Cm s KV( i - £ ^ ) 



The reversed voltage is now applied for a period t r * 

and considering the initial conditions that exist at 

time t 9 it follows that the equations for system 
s 

position and rate for t s Ik t - t^ are; 

{ 8 ) Ch) - VKA CM CM C m/< 

s*(s + A (s + \ /f) s(s * i/f) 5(5 + y#) 

and; 



( 9 ) C (*) s _ VK/tf C(tA 

5(5 + A (3 A) 



or; 



( 10 ) 



O.) •- VK 





and; 



(11 ) 



- VK 






For deadbeat response the system rate* C * 
must equal zero at time t- : , Therefore* by setting 
equation (11) to zero* the exact relationship between the 
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acceleration time,, t „ and deceleration time? t „ is 



obtained, 


namely? 


(12) 


.-/r _ £ - u « 

£ — /l- c 



To obtain the deceleration time required for any 
acceleration time, it is only necessary to take the 
natural logarithm of both sides of equation (12); 



(12a) 


i r * rJ>„(z-e' , '" f ) 



To obtain the system output for deadbeat response 
in terms of t g and t r , equation (12) is substituted 
into equation (10) which yields? 

( 13 ) CM ' VK(ti-ty) 

Finally, the system output in terms of the acceleration 
time only can be found by substituting equation (12a) 



Into (13) 


to give? 


(14) 


CM - VK\t> - tfJ.Z-C' /f \\ 


Thus 


by using a controller voltage proportional 



to R, and by choosing t s and t r as calculated from 
equations (12a) and (14), deadbeat response can be 
achieved for any step input provided that the plant 
Is adequately described by a (linear) second order 
transfer function, 

=> 1 8 = 
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4. A design example of the time invariant controller® 

A logical sequence of manipulations to determine 
the proper switching times for a given linear second 
order system can he easily obtained from equations 
(12a) and (14). For purposes of illustration*, assume 
that the system has the following characteristics? 

K = 10 Maximum output desired =10 

^=1.0 Saturation, will occur if KV is greater 
than 100® 

Using equation (14) with the constants given above 
and selecting a maximum applied controller voltage 
of 10 (to avoid the saturation non-linearity) a trans- 
cendental equation for the acceleration time*, t g9 is 
obtained? 

(14) iO = tOo[i* - &(z -(’**) 1J 

A value of t g equal to ®368 seconds will satisfy the 
equation. It is then necessary to find the deceleration 
time t r * To do this*, substitu** time t s into equation (12a) 

(12a.) i r ~ **( 2 ***• 

The total time of open loop operation is then ,6 36 seconds. 

Using the above constants*, the system equations were 
solved on the digital computer, Only one input was 
given to the system as other inputs would provide the 
same type of output changed by a constant only. Fig, 11 
shows the computed phase plane response of the system. 
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Test of the time invariant controller applied to 
a d-c motor plant* 

The controller analyzed in sections 3 and 4 above, 
is here Investigated experimentally in relation to a 
DC motor plant which was predominantly second order 
and linear up to the saturation limit of the driving 
amplifier • Runs were made with various step inputs in 
both the time” invariant and standard closed loop modes* 

The following operations by the open loop controller 
are required to test the theory: 

1) At time i c a step input is commanded* The 
loop is opened and an open loop driving voltage is 
applied to the plant* 

2) At time t s the driving voltage is reversed® 

3) At time the driving voltage is removed 
from the plant 9 the loop is closed 9 and the system 
returned to its standard closed loop mode of operation. 

The operations required were carried out by a timing 
device consisting of a variable speed D* C* motor 
geared to a shaft containing three “'make and break" 
switches. The time of contact of these switches rela~ 
tive to each other was variable® that is® their relative 
anglular positions on the shaft could be varied* The 
switches were connected to relay circuits which accomplished 
the desired switching* Thus® relay switching times 

could be varied by adjusting the speed of the D® C. 
motor or by adjusting the relative angular position of 
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the "make and break" contacts » 

Fig* 12 shows a block diagram of the experimental 
D. C* servo used for testing * When a step command is 
received 9 relay one switches to its normally open 
position* This is time t 0 » Relay two switches to its 
normally open position at time t g to reverse the open 
loop driving voltage* A third relay 9 not shown in the 
figure „ is used to return relays one and two to their 
normally closed positions at time t+« 




Fig, 12 Co Servo with a time invariant controller 
The transfer function of the open loop system was 
evaluated ass 

07 ) iJUJA 

$ (S + !/« 9) 

where the gain is 2 e 94 and the system time constant is 1.8 
seconds. The system output position^, C 9 was obtained by 
Integrating the output of a tachometer geared to the 
driving motor. As a result 9 the system output was 
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evaluated In volts g and consequently 9 the command and 
output signals could be compared In volts rather than 
radians* Further investigation of the open loop 
system revealed that an open loop command signal of 
about eight volts would cause saturation in the transistor 
power amplifier* A schematic of this amplifier is shown 
in Fig* 13® 

For the purpose of steady state accuracy 9 of 
primary importance in positioning sy stems 9 the closed 
loop gain was set at 29*4 or ten times greater- than the 
open loop gain* 

To simplify the construction of the controllers, it 
was decided to obtain switching times such that a one 
volt Input* R* would cause a one volt driving voltages, V. 
In other words* the constant of proportionality derived 
in equation (14) was made equal to one* This means that 
in equation (14)* V s C 9 and t g can be evaluated 
directly from that equation to yield the desired value* 



By using the value of t g equal to *98 sec* In equation 
(12a ) 9 t is found to be *54 sec® The relay switching 
times were then set accordingly* 

An input of three volts was used to check out the 



(14) 
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Fig. 13. Schematic of transistor amplifier. 
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switching time computed above* The response was not 
deadbeat* The swl tailing times were then adjusted to 
obtain deadbeat response to the three volt input* The 
actual time of voltage reversal was il 0 1 5 sec* compared 
to the computed *98 see* The total time of operation 
was 1*55 sec* compared to the computed time of 1 e 52 sec* 
These differences were probably due primarily to the fact 
that the system was only quasi linear and second. order 8 
and the 4 #^**** j Lant were determined 

assuming a linear second order systems 

A six channel Brush recorder was used to record 
the closed loop command signal 8 the ©pen loop command 
signal^, the system output and the system rate* Runs 
were made with the inputs varying from 1*4 volts to 10 
volts* A standard closed loop response to the same 
step commands was also made for comparison purposes* 
Figures 14 through 21 show the responses obtained* Each 
page contains a tlme° invariant system response followed 
immediately by a standard closed loop response t© the 
same signal* Table I is a comparison of the different 



runs 
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!o M Km* ,-. i maximum output to desired output* 

pi 

2 o Time »c ate a ay state operation, , taken to be 
when ?*y .•?•*£!» i« within *t volts' of command 
signal « 

3® ^ tg<, Is defined as the difference between 
the tira«p *c rte-ady state for the epen=> 
closed I'. ,p and ^1" sed lo;*p operation* 



Deadbeat reepons^e or near it^i' fc*jat response;, 
could be force - ng th* *. invariant controller 
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Fig* l4. Time invariant and closed loop response of DC-servo 

system with a 1*94 volt input* - 
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Pig* 15* Time invariant and closed loop response of DO servo 

system with a 5*^5 volt input* 
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'Pig* 16. Time invariant and closed lo-‘P response of DC servo 
system with a 4.0 volt input* 
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Fig. 17* Time invariant and closed loop response of DCLeervo. 
system with a 5* 0 volt input® 




Fi g* 18* Time invariant and closed loop response -of LC servo — 

system with a 6*p volt input * 
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Fig* .21* Time invar iant^and closed loop, response of DO servo., 
system with a 1.4 volt input* 
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The controller reversed the Input to the motor at the 
proper time but the reversal Is not indicated on this 
trace. It was not possible to record the Input to 
the motor (amplifier output) showing the reversal of 
driving voltage,, because of grounding problems,, A 
ground at the motor input would cause one or both of 
the power transistors to burn out in the power ampli- 
fier. The third trace is the system output and the 
fourth trace, the system rate. The time t s , the time 
at which the polarity or sign of V is reversed^ is 
readily observed on the last trace, system rate 9 as 
there is a definite change in slope at that time® 

Time on all the traces reads from left to right with t y 
signified by the commencement of the open loop command 
signal and t t by the end of this command and the start 
of the closed loop command signal. 

It was mentioned earlier that an output signal of 
eight volts would saturate the power amplifier It 
follows, then, that if , the amplifier will 

saturate when the driving voltage is reversed. The 
reason for this is that at the time of voltage reversal 
a back electromotive force has built up in the motor 
which approaches the driving voltage. This back emf 
has the effect of increasing current requirement from 
the amplifier when the driving voltage is reversed. 
Therefore, at the time of switching, the system should 
tend toward saturation if the signal voltage is greater 
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than four volts* This saturation caused a reduction 
in the braking or deceleration power available* Prom 
this it follows that the system output would tend to 
overshoot. These effects may be seen in Figures 16, 

17 , and 18. By analyzing the system rates at time t s 
and t.j. for these figures it can be demonstrated that 
saturation does occur in the time increment t s S t ‘s t t « 

At time t s , for the three figures mentioned above, the 
system rate increased proportionately with the command 
signal as called for by the theory. Yet at time 
the system rate was larger than predicted meaning that 
the system did not decelerate sufficiently which 
resulted in overshoot at time t^. For example, the 
system rates at time t g and for the five volt input 
(Fig. 17) are 6.1 and .95 volts per second respectively. 
For an input of 6.2 volts (Fig. 18) one would predict 
by the principle of superposition that the rates would 
be increased by the factor 6. 2/5.0 or would be 7®5 
and 1.1 volts per second. Actually, Fig. 18 shows 
the rates to be 7.5 and 1.5 volts per second. The 
rate at time t g agrees with the theory but at time t^ 
the rate differs by .4 volts per second. This difference 
can be accounted for by the saturation that occurs 
when the voltage Is reversed. 

When the system Input exceeds eight volts, saturation 
is predicted at t 0 and at t s . One would expect, the re = 
fore, that the system output would never reach the 
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desired point at time t^. or in other words s the system 
will have undershoot at time t t » The reason for this 
phenomenon is that the driving voltage Is never large 
enough to provide the response called for by the 
theoretical equations. The results of this type of 
saturation are shown in Figures 19 and 20. Fig. 20 
clearly demonstrates the effect on system output. 

Fig. 21 shows the response of the system with a 
small input, 1.4 volts. Here the system output*, by 
time t t , does not reach the desired value. The primary 
cause for this is that the small driving voltage is of 
insufficient magnitude to overcome the stietion and 
friction in the system. 

From these experimental tests it can fee seen that 
the time invariant switching scheme does provide 
deadbeat response over a range of inputs. If deadbeat 
response is desired in a positioning feedback control 
system, the method of time invariant switching can 
provide this in a given time increment over a designed 
range of inputs. 



6. Test of the time invariant controller applied to 

an amplidyne driven motor. 

The time invariant controller has been shown to 
work over a range of Inputs with a quasi linear , 
second order system. One might now inquire as to its 
usefulness with a system that is quite nonlinear in that 
it exhibits hysteresis, excessive stietlon, and satura- 
tion. To investigate this, a system was constructed 
that incorporated an amplidyne, a i HP d“C shunt wound 
motor, and a permanent magnet, D, C. generator. The 
output of the generator, or tachometer, was integrated 
to simulate the system output. As in the previous 
experiment the system output and input was in volts. 

Pig, 22 depicts the system as set up in the laboratory 
with the relays performing the same functions as out- 
lined in section 5» Pig, 23 is a picture of the system 
as it appeared in the laboratory. 




Pig. 22 Amplidyne driven motor servo system with 
time Invariant controller. 
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A plot of the voltage Input versus the tachometer 
voltage output for the open loop system appears on 
Pig. 24. Notice that at least two volts are required 
to start the system in motion and that at about six 
volts the system starts to go into saturation. The 
hysteresis loop is rather broad which could provide 
difficulty in system performance. Prom open and closed 
loop tests of the system* without the time invariant 
switching* the open loop transfer function was determined 
to be approximately? 



(17) 

S(S t 7) 

Depending on the test being conducted* the system gain 
varied from 100 to 115 and the reciprocal time constant 
from 6.9 to 7.5 inverse seconds. 

Using the values of gain and the time constant 
from equation $7) and having a unity relationship between 
the open loop output and open loop driving voltage* t_ 
was determined from equation (14) to be about ,33 
seconds. Then t^ g from equation (12a)* was evaluated 
as .27 seconds. The system was rim with an input step 
of four volts. The response was not deadbeat so the 
times of switching were adjusted to yield a near dead= 
beat response. The timing circuit controls were not 
fine enough to permit precise adjustments of the switching 
times. The adjusted times were then measured and found 
to be .38 second for time t g and .09 second for t r . 
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Fig. 24. Hysteresis loop of aaplidyne driven servo systen. 
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The acceleration time was .05 second longer than the 
computed value and the deceleration time was .18 second 
shorter than the computed value. This shorter decelera- 
tion time is probably due to the large friction that 
was present in the motor^ 

A Brush recorder was used to record the values of 
the closed loop command signal, the open loop command 
signal, the system output, the system rate, and the 
voltage applied to the motor. These recordings are 
shown in Figures 25 through 38* Following each time® 
invariant run to a given input, the system was run 
in the normal closed loop mode with the same quantities 
being recorded. A figure of the closed loop recording 
follows each time invariant figure. The maximum over® 
shoot and time~to®steady®state was evaluated for all 
runs and a comparison of the times to steady state was 
made. The results of this evaluation are tabulated in 
Table II. 
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£*«*tt*«««**«*####** 

Amplldyne * Open-closed loop * Closed loop 

* * S l r J°» 

* Input * M pt ' tgg 3 - * M pt t gB At 3 



* 



3.0 volts 


1 .0 


1.18 sec. 


1 .0 


.59 sec. 


<=.59 se 


3.5 


1 oO 


1 .00 


1.0 


.50 


°.50 


4.0 

* 


1 .0 


.50 


1.26 


1.00 


.50 


4.5 

* 


1 .0 


CO 

o 


1.31 


1.05 


.57 


5.0 


1 ,0 


.48 


1.58 


1 .00 


.52 


5.5 


1 .25 


1 .03 


2.02 


.93 


<=. 1 0 


6.0 


1 .88 


1 .00 


2.12 


1.30 


.3 




* # 4* # * 


* * # # it # 




# # * # # 


Vt 



Note: 

<• v : Ratio of maximum output to desired output. 

2. Time to steady state operation taken to be 
when system is within .1 volts of command 
signal. 



3. 



At gs is defined as the difference between 
the times to steady state for the open- closed 
loop and closed loop operation. 



Prom the recordings it can be seen that the range 



from four to five volts the response was almost dead- 
beat. (Figures 25 through 30.) When the input voltage 



was in excess of five volts the system output started 



to overshoot the desired value at time t^. 
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Fig. 2 6* Closed loop response of amplidyne driven motor servo 
system with a four volt input. 
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Fig* 27* Time invariant response of amplidyne driven motor servo 
_ system with a volt input. 
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Fig. 2<3. Closed loop response of amplidyne driven motor servo 
system with a 4-jy volt input. - - 
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Fig. 29* Time invariant response of amplidyne driven motor servo 
system with a five volt input.- 
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Fig. 50. Closed loop response of anplidyne driven motor servo 
system- with a five volt— input*. 
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section 5 with the D. 0 o servo 9 this is due to satura® 
tion of the amplidyne when the driving voltage is 
reversed. This type of response is shown in Figures 31 
through 34. 

For step inputs less than four volts system 
performance deteriorated rapidly. This was due almost 
entirely to the excessive stiction and friction in the 
system. It was found that the system did not respond to 
a command less than two volts or about 30 % of the 
operating range. The system response to inputs of 3 
and 3.5 volts are shown in Figures 35 through 38 . 

As might be expected » the use of an open^loop 
controller with a far-from= ideal plant does not lead 
to acceptable performance. 
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Fig. 51. Time invariant response of amplidyne driven motor servo 
system with a 5-§ vdlt input. 
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Fig. 52. Closed loop response of amplidyne driven motor servo 
system with -a- 5a' volt input. 
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Fig* 5> Time invariant response of amplidyns driven servo 
system with a six volt input* 
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Fig. j4. Closed loop response of amplidyne driven servo 
_ system with a six volt input. 
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Pig. 55* Time invariant response 
system with three volts input. 
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Fig* 56 * Closed loop response of amplidyne driven motor servo 
L : system-with-a -three volt- input* 
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Fig. 57. Time invariant response to amplidyr.e driven motor servo 
system with a 5 e‘ volt input. 
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Fig. 53. Closed loop response of amplidyna driven motor servo 
system with a ^ volt input. 
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7. Conclusions. 

In the search for a simple scheme to provide 
deadbeat response for step inputs to a second order 
system, attention was concentrated on an open-closed 
loop controller. For the system investigated, switching 
time in the open loop mode were fixed and the output 
from the controller to the plant was proportional to 
the magnitude of the step input. 

Using this scheme, any linear system can be 
controlled with deadbeat response (Principle of super- 
position) . To control the second order plant considered 
in this study only one switching point is required. 

The experimental program was intended to suggest the 
shortcomings of this approach when applied to a relay 
system with its non-linearities and higher order dynamics. 

The first real system studied was relatively free 

from any non-linearities and provided near deadbeat 

response over approximately 50$ of the designed operating 

range. A system with rather pronounced non-linearities 

was also chosen for study and provided near deadbeat 

response over less than 20$ of the desired operating 

range. It is concluded, then, that if deadbeat or near 

deadbeat response is desired, the time invariant controller 

investigated in this thesis can be used to force this 

response in a linear system . With a non-linear system, 

the controller is of little use. In any event, the time 

of response to any size step requires the same Interval 

of time for open loop operation. 
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APPENDIX I 



PROGRAM ANALOG 

A Control Data Corporation CDC 1604 digital 
computer is available at the U* S. Naval Postgraduate 
School. Dr. J. R« Ward of the Postgraduate Schools 
Department of Electrical Engineering 9 has written a 

O 

FORTRAN language program for this computer which 
simulates an analog computer. The program is referred 
to as Program Analog* It retains all the versatility 
of the analog computer and provides both a graphical 
and numerical output. With this program many of the 
problems of magnitude and time scaling as well as the 
difficulties in simulating nonlinearities associated 
with an analog computer are avoided. The program „ which 
is reproduced in Figures 39 through 41 8 uses a modified 
Runge-Kutta integration scheme* 

This program was used for the investigations 
carried out in section 2* Fig. 39 is the main part of 
the program. It consists of the necessary Instructions 
for the analog simulation of a given problem and does 
not require any changes as the problem is varied. This 
part of the program may be thought of as the analog 
computer. Fig. 40 is the FORTRAN simulation of the 

^D. D, McCrackens, A guide in FORTRAN programmings, 
John Wiley & Sons*, Inc. 8 New York 9 London^ 1961. 
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problem board# This part of the program, must be 
changed in the same manner as one would change the 
wiring on a problem board* Fig* 41 » the data input 9 
sets the initial conditions^ potentiometer settings^ 
problem time and provided means by which various output 
data can be obtained* A typical numerical output 
for a Program Analog problem appears in Pig* 42* 

Figures 3 through 10 in the body of the these o re 
typical outputs for the program* 
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.JOB MALLEY ( 10 MINUTES SHOULD CO) 

PROGRAM ANALOG 
CALL EXEC 
ENO 

SUBROUTINE EXEC 

DIMENSION X(15C),XCOT(100),C(15C),ITITLE(10),JTITLE(10),KTITLE(40) 

1 ,IT(5),IP(10),IG(10),CT(5),XO(100),PR(10),GR(1C),X1 (900 ) , 

2 Y1 (9 0C ) ,X2 ( 900 ),Y2(900),X3(9C0),Y3(900),X4(9C0),Y4(900), 

3 X5 (90C ) ,Y5 ( 900 ) ,A ! 100 ) 

DO 1 J = 1 , 1 50 

X ( J ) = 0. 

1 C ( J ) = 0. 

DO 1C01 J=1,10C 

1001 X0( J ) = 0. 

DO 1002 J = l, 1C 
I T I TLE ( J ) = 8H 

1002 JTITLE(J) = EH 
DO 1003 J= 1 ,40 

1003 KT I TLE ( J ) = 8H 
' NRC = 0 

PRINT 200 
200 FORMAT ( 1H1 ) 

2 C 1 FORM ATI / ) 

202 FORMAT!//) 

203 FORMAT! ///) 

READ 102, ( I T ITLE! J ) , J = 1 ,5) 

102 FORMAT! 10A8) 

READ 101, N, ITEST 
101 FORMAT ( 12, A8 ) 

ICHECK = 8 H EOUATIC 
IF! ICHECK - ITEST) 2,3,2 

2 PRINT 204 

204 FORMAT! 18H DATA FORMAT ERROR) 

STOP 1 

3 READ 101 , NR, ITEST 
ICHECK = 8H RUNS 

IF! ICHECK - ITEST )4, 6, 4 

4 ICHECK = 8H RUN 

I Ft ICHECK - ITEST) 7,6,7 

6 I F ( N R - 9) 3,8,7 

7 PRINT 204 
STOP 2 

8 READ 101 , NT, ITEST 
ICHECK = 8H TITLE C 

IF! ICHECK - ITEST) 10,9,10 

9 I F ( NT - 9)11,11,10 

10 PRINT 204 
STOP 3 

11 DO 12 1=1, NT 

READ 102, ( JTITLE(J),J=1, 10) 

12 PRINT 205, ( JTITLE! J),J=1 , 1C) 

I F ( N R - 1)10,13, 14 

13 PRINT 206 
GO TC 15 

14 PRINT 207, NR 

205 FORMAT ( 5 X , 1CA3) 

206 FORMAT! 5X.20F1 RUN IS CALLED FOR.) 

207 FORMAT (5X, II ,21F RUNS ARE CALLED FOR.) 

15 PRINT 203 
300 NRC = NRC + 1 

READ 102, ITEST 

I CHECK 1 = 8HZFRC CCE 

ICHECK2 = 8HFCLC CCE 

IF! I CHECK 1 - ITEST) 16, 18 , 16 

16 I F ( I CHECK2 - I TE ST ) 1 7 , 20 , 1 7 

17 PRINT 204 
STOP 4 

18 DO 19 J=1 , IOC 

19 C(J) = 0. 

20 RE AC 10 1 ,NC, ITEST 
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21 , 22,21 



, K= 1 , 5 ) 



ITEST) 52,32,32 



ICHECK = 3H CTEFF. 

I F ( ICHECK - ITEST) 

21 PRIM 2C4 

STOP 5 . 

22 I F (NC) 23, 2d, cZ 

23 DO 24 J=1,MC 
RE AC 103, ( IT ( K ) » CT ( K ) 

103 FORMAT (5(14, F12. 4)) 

DC 24 Ks 1 , 5 
ITK = I T ( K ) 

24 C( ITK) = CT ( K ) 

25 PRINT 203, NRC 

208 FORMAT (12H RLN NUMBcR 
IRE AS FOLLOWS,/) 

K = C 

DC 27 J=l,99 
I F (C ( J ) )26, 27 , ct 

26 PRINT 2C9, J , C ( J ) 

209 FORMAT ( 1CX2HC (, 12, 4H ) = 

K = K + 1 

27 CONTINUE „ 

IF(K)1027, 10^6,1 027 

1026 PRINT 240 

240 FORMAT ( 5X4HNCNE ) 

1027 PRINT 2C2 

READ 102, ITEST „ r r c 
I CHECK 1 = 8PZt_RC ICS 
ICHECK2 = 3HFCLC ICS 
ICHECK3 = 8HFCLC FCS 
I F ( I CHECK 1 - IT EST ) AC , 28 , 2C 
23 DC 29 J = 1 , N 

29 XC(J> = 0 . 

30 F ( I CHECK 2 - ITEST )3 1^35 ,31 ^ 

31 I F ( ICHECK3 - 

32 PRINT 204 
STOP 6 

33 DO 3 4 J = 1 » N 

34 XC ( J ) = X ( J ) 

35 READ 101 , N I 
ICHECK = 3H 
I F ( I CHECK - 

'36 PRINT 2C4 
STOP 7 

37 I F t N 1)38, 40,33 

38 DO 39 J = 1 ,NI 
READ 103, ( IT ( K ) 

DC 39 K = 1 , 5 
ITK = I T ( K ) 

39 XC( ITK) = CT (K) 

40 PRINT 2 1 C 

210 FORMAT (48H THE 
K = C 

DO 42 J=1,N 
IF ( XC( J ) ) 41 , 42 , 4 1 

41 PRINT 211 , J,XC( J 

211 FORMAT ( 1 0X3HXC ( , I 2 , 4H) - 

K = K + 1 

42 CCNT INUE 
1 = ( K ) 1042 , 1 041 , 1 042 

1041 PRINT 240 
1 C42 PRINT 202 

READ 1C2,iT:_T 
I CHECK 1 = BHHCLC TI v 

ICHECK2 = 8HPEAC TI-- 
I F ( I CHECK 1 - ITcST 43,46 

43 I F I I CHECK2 - ITcST)44, 

44 PRINT 204 
STOP 8 

45 READ 104,TC,Tr 



I 1 , / / / , 46H THE NON-ZERC DATA COEFFICIENTS A 



E12.5) 



ITEST 
IC CARC 
ITEST ) 36 



37, 36 



, CT (K ) , K= 1 , 5 ) 



NCN-ZERO INITITAL CONDITIONS ARE AS FOLLOWS,/) 



E12.5) 



43 
.5 ,44 
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1 04 FORMAT ( 10F8. 4 ) 

46 RLAD102 , ITEST 
ICHECK! = 8HCCMPUTE 
ICHECK2 = 8HF0LC STE 
ICHECK3 = 8HREAC STE 

I F ( I CHECK 1 - ITEST )43, 47 , 48 

47 INDIC1 = 1 

READ 101, KP, ITEST 

I CHECK = 8H IS CRDE 

I F ( I CHECK - ITEST ) 50 , 59, 50 - 

48 I F ( ICHECK2 - I T E ST ) 4 9 , 59 , 4 9 

49 I F ( I CHECK3 - I TE ST ) 5 0 , 51 , 50 

50 PRINT 204 
STOP 9 

51 INDICT = 0 

READ 104 , (C ! J) , J = 1 C3 , 109 ) 

I F ( C ( 10 3) - TO 152,53,52 

52 PRINT 204 
STOP 10 

53 I F ( C ( 104) 154,52,54 

54 NOT = 1 

I F ( C ( 1 0 5 1 - TF 1 55 , 59 , 55 

55 I F ( C (106) 156,52,56 

56 NDT = 2 

I F ( C ( 1 07 1 -TF)57,59,S7 

57 I F ( C ( 1 0 8 1 158,52,58 

58 NDT = 3 

I F ( C ( 109) -TF) 52, 59, 52 

59 PRINT 212, TO , TF 

212 FORMAT { 3 1 H TPE TIMING DATA ARE AS FOLLOWS ,// 

1 5X, 15FINI TIAL TIME = ,E11.5,/ 

2 5X, 151-FINAL TIME = , E 1 1.5) 

IFIINDIC1 - 1)61,60,61 

60 PRINT 213, KP 

213 FORMAT { 5X, 35HTHE STEP SIZE IS COMPUTED, BASED ON ,/ 

1 5X , 31 FSMALLEST VARIABLE CF CRCER l.OE ,12) 

AP = 10.0**KP* 1 .0E-04/ITF - TO 
DO 1060 J= 1 , N 
1060 A ( J ) = AP 

A ( 1 00) = TF -TC 
C { 1 04 ) = All CO 1*1. CE-05 
GO TO 63 

61 NDTT = NDT«2 

DO 62 K = 1 , NDTT, 2 

62 PRINT 214, C( 1C3+K) ,C ( 102+K) ,C( 104+K) 

214 FORMAT! 5X, 15HSTEP SIZE = ,E11.5,10H FROM T = ,E11.5,3H TO T 

1 El 1.5) 

63 PRINT 202 
READ 102, ITEST 
ICHECK = 8HHCLD PRI 

I F ( I CHECK - ITEST) 64,74, 64 

64 ICHECK = 8HREAD PRI 

IF! ICHECK - ITEST)65,1064,65 
1 064 READ 101, NP, ITEST 

ICHECK = 8H VARIABL 
IFfICHECK- ITEST ) 65 ,66, 65 

65 PRINT 204 
STOP 11 

66 IF(NP)67,74, 67 

67 I F { N P - 10)68,68,65 

68 READ 1104 { I F t K ) ,K = 1 ,NP) 

1104 FORMAT (1014) 

READ 101, NT, ITEST 
ICHECK = 8H TITLE C 
IFfICHECK - ITEST)65,69, 65 

69 IF(NT)70, 73,70 

70 I F ( N T - 4)71 ,71 ,65 

71 DO 72 1=1, NT 

K = (I - 1)*10 +1 
KP9 = K + 9 
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72 READ 102, (KTITltl J),.J = K,KP9) 

73 RE AC 10 1, INC P R , I TEST 
ICHECK = 8H INCREME 
IFIICHECK - ITEST >65,74, 65 

74 READ 102, ITEST 
ICHECK - 8FHCLD GRA 

I F { ICHECK - ITEST)75,81, 75 

75 ICHECK = 8HREAD GRA 

IF ( ICHECK - ITEST) 76 ,77, 76 

76 PRINT 204 
STOP 12 

77 READ 101 ,NG, ITEST 
ICHECK = 8H GRAPHS 
IFIICHECK - ITEST ) 76 ,78, 76 

78 IFING-) 79, 81*79 

79 IFING - 5)80,80,76 

80 NG2 = NG*2 

READ 1 104, { IGIK ) ,K=1 ,NG2 ) 

READ 101 , INCGR, I TEST 
ICHECK = 8H INCREME 
IFIICHECK - ITEST ) 76 ,8 1, 76 

81 PRINT 215 

215 FORMAT! 17H PRINT SUMMARY ,/) 

IF(NP)83,82,S3 

82 PRINT 216 

216 FORMAT ( 5X , 1 1 FNO PRINTOUT) 

GO TO 84 

83 PRINT 217, INCPR 

217 FORMAT (5X, I 2 ,29h INCREMENTS BETWEEN PRINTOUTS, 

1 5X , 25FTHE VARIABLES PRINTEC ARE,/ ) 

PRINT 218, I IPI J ) ,J = 1 ,NP) 

218 FORMAT! 10X,2FX{ , 13, 1H) ) 

DO 1 C84 J = 1 , N P 

IF! IPI J) ) 1084, 1 C 85 , 1084 

1084 CONTINUE 
GO TO 84 

1085 PRINT 1217 

1217 FORMAT !/,5X,22H* X ( 0 ) REPRESENTS TIME ) 

84 PRINT 202 
PRINT 219 

219 FORMAT! 17H GRAPH SUMMARY ,/) 

IF(NG)86,85, 66 

85 PRINT 220 

220 FORMAT I 5X, 9FN0 GRAPHS) 

GO TO 2089 

86 PRINT 221 , INCGR 

221 FORMAT! 5X, 12, 26H INCREMENTS BETWEEN POINTS,/ ) 
PRINT 241, I G ! 1 ) , I G ( 2 ) 

241 FORMAT ( 1 OX , 1 3HGRAPH A IS Xt,I3,8H) VS. X ( , I 3 , 
IF! NG - 2)87,1087, 1087 

1087 PRINT 242, IG I 3 ) , IG(4) 

242 FORMAT ( 1 OX , 1 3HGRAPH B IS Xl,I3,8H) VS. X { , I 3 , 
I F ! NG - 3)87,1088, 1088 

1088 PRINT 243, I G ! 5 ) , IC!6) 

243 FORMAT { 1 OX, 13HGRAPH C IS X!,I3,8H) VS. X ( , I 3 , 
I F ( NG - 4)87, 1089, 1089 

1089 PRINT 244, I G ( 7 ) , I G ( 8 ) 

244 FORMAT ( 1 OX , 1 3HGRAPH C IS X(,I3,8H) VS. X!,I3, 
I F ( N G - 5)87,1090,1090 

1090 PRINT 245, IGI9) , IG! 10) 

245 FORMAT ( 1 OX, 1 3HGRAPH E IS X!,I3,8H) VS. X ! ,13, 

87 DO 2087 J=1 ,NG2 

I F ( I G ( J ) ) 2 08 7, 2088, 2087 

2087 CONTINUE 
GO TO 2089 

2088 PRINT 1217 

2089 PRINT 200 

PRINT 205, ( ITITLE! J), J = 1 ,5) 

PRINT 223, NRC 

223 FORMAT (5X, 10FRUN NUMBER, I 2 , / // ) 



1H) ) 
1 H ) ) 
1H) ) 
1H ) ) 
1H) ) 
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IF(NP)88,90, 88 

88 I F (NT)89,90, 89 

89 PRINT 224, (KTITLEf J) ,J=1 ,KP9 ) 

224 FORMAT ( 2 X , 9 ( A 8 , 4X) ,A3) 

PRINT 201 

90 T = TO 

DT = C( 104 ) 

DO 91 J=1 ,N 

91 X ( J ) * XO ( J ) 

LINES = 0 
NCPTS = 0 
NUMPTS = 0 

301 IF(NP)302, 312,302 

302 IF(N0PTS)303,308,303 

303 I F ( XMCDF(N0PTS»50*INCPR) 1304,306,304 

304 I F ( XMCDF ( NOP TS , 10* INC PR) 130 5,307,30 5 

305 I F ( XMODF ( NOPTS t INCPR) 1312,308,312 

306 PRINT 200 

307 PRINT 201 

308 CALL DERIVI T,X, XDOT,C) 

LINES = LINES + 1 

00 311 J = 1 ,NP 

1 F ( I P ( J 1 1310 ,309,310 

309 P R ( J 1 = T 
GO TO 311 

310 IPJ = I P ( J 1 

P R ( J 1 = XI IPJ) 

311 CONTINUE 

PRINT 225, ( PR ( J ) , J= 1 ,NP) 

225 FORMAT ( 1 0 ( IX, El 1.5)1 

312 IFING1313, 31 8,313 

313 I F ( XMCDF (NOP TS, I NC GR ) ) 3 1 8 , 3 1 4 . 3 1 8 

314 I F( XMODF (NOPTS, IMCPR) 1 1315, 1314, 1315 
1315 CALL DERIV(T,X,XDOT,C) 

1314 DO 317 J=1 , NG2 

I F ( I G ( J 1 1316,315,316 

315 GR ( J 1 = T 
GO TO 317 

316 IGJ = I G ( J 1 

GR ( J 1 = X( IGJ) 

317 CONTINUE 

NUMPTS = NUMPTS + 1 
Y 1 (NUMPTS) = GR ( 1 1 
XI (NUMPTS) = GR( 2) 

Y2( NUMPTS) = GR ( 3 ) 

X2 ( N UMP TS ) = GR ( 4 ) 

Y 3 ( NUMPTS) = GR ( 5 ) 

X 3 ( NUMPTS) = GR ( 6 ) 

Y4 (NUMPTS) = G R ( 7 ) 

X4 ( NUMP T S ) = GR ( 8 ) 

Y5( NUMPTS) = GR ( 9 ) 

X5 ( NUMP TS ) = GR( 10) 

318 NOPTS = NOPTS + 1 

I F ( L I N E S - 250)1319, 1318,131 8 

1318 PRINT 1216 

1216 FORMAT (//24F STOP AT 250 PRINT LINES) 
GO TO 341 

1319 I F ( T - 1. £ + 04)320, 319, 31 9 

319 PRINT 226 

226 FORMAT ( / / 1 9H STCP AT T = 10,000) 

GO TO 341 

320 IFtNOPTS - 1C000 )322,321,321 

321 PRINT 227 

227 FORMAT ( //26H STCP AT 10,000 INCREMENTS) 
GO TO 341 

322 I F ( NUMPTS - 900)324,323,223 

323 PRINT 228 

228 FORMAT ( //25H STCP AT 900 GRAPH POINTS) 
GO TO 341 

324 I F ( T - TF) 326,325, 325 
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325 PRINT 229 

229 FORMAT ( / /26H NORMAL STOP AT FINAL TIME) 

GO TO 341 

326 DO 328 J=1 , 1 50 

I F ( A8SF (X ( J) ) - 1. £ + 04)328,327, 327 

327 PRINT 230, J 

230 F0RMAT(//1 lH STOP AT X(,I3,1CH) = 1C, COO) 

GO TO 341 

328 CONTINUE 

330 I F ( I NO I C 1 >331. 332* 1 331 

331 STOP 13 

332 CALL RKUTTA (N,T,DT,X,C) 

I F ( NOT - 1 )334, 333,334 

333 DT = C ( 1 04 ) 

GO TO 340 

334 I F ( NOT - 2) 337,335,337 

335 I F ( T - C ( 105) )333, 336,336 

336 DT = C ( 106 ) 

GO TO 340 

337 I F ( T - Cl 105) 1333,338,338 

338 I F ( T - C( 107) )35o,339,339 

339 DT = C ( 108) 

340 T = T+DT 
GO TC 301 

1331 CALL RKUTTA2 ( N , T , CT , X , A , C ) 

GO TO 301 

341 I F(NG) 1341,356,1341 

1341 GO TO ( 342,343,344,345,346,347,348,349,350) ,NRC 

342 ITITLE16) = £H RUN 1 
GO TO 351 

343 ITITLE16) = 8H RUN 2 
GO TO 351 

344 ITITLE16) = 8H RUN 3 
GO TC 351 

345 ITITLE16) = 8H RUN 4 
GO TC 351 

346 ITITLEl 6) = £H RUN 5 
GO TO 351 

347 ITITLE(o) = 8H RUN 6 
GO TO 351 

348 ITITLE16) = £H RUN 7 
GO TO 351 

349 ITITLE16) = 8H RUN 8 
GO TO 351 

350 ITITLEl 6 ) = 8H RUN 9 

- 35 1 LABEL = 4H 

MCDCLRV = 0 
SFX = 0. 

SFY = 0. 

MINCFFX = 0 
MINOFFY = 0 
LABELNO = 11 
MODE = 0 

ITITLEl 7 ) = 8H GRAPH A 

CALL GRAPH2 ( NUN PT S , X 1 , Y 1 , 8 , MCDCURV, L ABEL , I T I T L E , SFX , SF Y, 
1 MINCFFX, MIN0FFY,LA8ELNC»MCDE) 

I F ( N G - 2)356,352,352 

352 ITITLE17) = EH GRAPH B 

CALL GRAPH2 { NUM P T S , X2 , Y 2 , 8 , MODCURV, L ABEL , I T I TLE , SFX , SF Y, 
1 MINOFFX, MINOFFY, LABEL NC, MODE) 

I F ( NG - 3)356,353,355 

353 ITITLEl 7) = EH GRAPH G 

CALL GRAPH2 ( NUM PT S , X 3 , Y 3 , 8 , MCDCURV , L ABEL , I T IT Lc , SF X , S F Y, 
1 MINOFFX, MINCFFY,LABELNC,MCDE) 

IF(NG - 4) 356,354,354 

354 ITITLE17) = cH GRAPH D 

CALL GRAPH2 ( NUMPT S, X4 , Y4 , 8 , MCDCURV, L ABEL , I T I TLE , SFX, SF Y, 
1 MINOFFX, MIN CFFY,LABELNC, MODE) 

I F ( NG - 5)356,355,355 

355 I T I T L E ( 7 ) = EH GRAPH E 
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CALL GRAPH2 (NUMPTS, X5 , Y 5 , 8 , MODC URV , LABEL, [TITLE, SFX, SPY, 

1 MINQFFX,MINCFFY , LABELNC , MODE ) 

356 PRINT 200 
I F ( NRC - NR) 300,357, 357 

357 RETURN 
END 

, SUBROUTINE RKUTTA (N,T,DT,X,C) 

DIMENSION X( 150) , AK{4, IOC) ,X COT (100) , XC( 100) ,CT( 4 ) , C (150 
CT ( 1 ) = 0.0 

C T ( 2 ) = 0.5 
CT ( 3 ) = 0.5 
C T { 4 ) = 1.0 
DO 4 1=1,4 
TC = T + CT ( I ) *0T 
DC 2 J= 1 , N 

2 XC(J) = X ( J ) + CT( I)*AK( 1-1, J) 

CALL DERIV (TC,XC,XDOT,C) 

DO 4 J= 1 , N 

4 AK ( I , J ) = DT «XDCT ( J ) 

DO 3 J = 1 , N 

3 X(J) = X(J) + (AK( 1,J)+2.*AK(2, J)+2.*AK(3,J) +AK(4,J ) )*0.1 6666 66667 
RETURN 
END 

SUBROUTINE RKUTTA2 ( N , T, CT , X , A , C ) 

THIS SUBROUTINE IS BASEC ON THE METHOD DESCRIBED IN COMM. OF 
THE A.C.M., JUNE I960 (PP. 355 -360). IT CHECKS TFE 
TRUNCATION ERROR AT EACH STEP OF THE INTEGRATION AND ACJUSTS 
THE STEP SIZE ACCORDINGLY. THE ACTUAL RUNGE-KUTTA INTEGRATION 
IS PERFORMED EY SUBROUTINE RKUTTA, WHICH IS AVAILABLE IN THE 
USNPGS COMPUTING CENTER. THE ARGUEMENTS ARE, 

N = NUMBER OF (FIRST ORCER) EQUATIONS. MAXIMUM N =99. 

T = TIME AT START OF INTEGRATION STEP (UPDATED BY THIS 
SUBROUTINE AFTER THE COMPLETION OF EACH STEP). 

DT = STEP SIZE AT START OF EACH STEP (ALSO UPDATED). 

X ( I ) = THE N DEPENDENT VARIABLES (ALSO UPDATED). 

A(I) = THE SPECIFIED ALLOWABLE ERROR PER UNIT TIME FOR EACH 
OF THE DEPENDENT VARIABLES. 

A ( 100) IS USED TC ENTER THE TOTAL TIME. 

NOTE THAT IF NECESSARY DT IS REDUCED FROM THE VALUE STATED 
IN THE ARGLEMENT UNTIL THE SPECIFIED ACCURACY HAS BEEN 
ACHI EVEC. 

DIMENSION X ( 150 ) ,XS( ICO) , X22 ( 1 CO ) , X2 ( 1 00 ) , A ( 1 C 0 ) , C ( 1 50 ) 

TS = T 
H = 2.0*DT 
DO 1 I = 1,N 
X2( I ) = X ( I ) 

1 X22 ( I ) = X ( I ) 

CALL RKUTTA ( N , T S , H , X2 2 , C ) 

H = DT 

2 CALL RKUTTA (N,TS,H,X2»C) 

DO 3 1=1 , N 

3 XS ( I ) = X2 ( I ) 

TS = TS + H 

CALL RKUTTA ( N , T S , H , X S ,C ) 

U2 = 0.03 
DO 6 1=1, N 

E 2 1 = { X22 ( I ) — X S ( I ) )*0. 06666666667 

E 2 1 R = DIMF ( ABSF ( E21 ) , AB SF ( ( ABSF ( XS ( I ) ) +ABSF ( X22 ( I ) ) 

1 -1 .99*ABSF (X ( I ) ) ) *1 .0E-C9 ) ) 

C THIS CONDITION PREVENTS RCUNC-OFF ERROR FROM TAKING CONTROL. 

U = E21R / ( A ( I)*6.0*H) 

IF (U - U2 ) 6 , 6 , 5 

5 U2 = U 

6 X S ( I ) = XS ( I ) - E21 
IF ( U2 - 1.0)11 ,7,7 

7 IF (H - A ( IOC )* 1 .0E-09 )8 ,8 ,9 

8 DT = A( 100)*1.0E-09 

C THIS SETS THE MINIMUM STEP SIZE TO 1.0E-09 TIMES THE TCTAL TIME. 

GO TC 16 

9 DO 10 I = 1 ,N 
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X22 ( I) = X 2 { I ) 

10 X 2 t I ) = X ( I ) 

TS = T 

H = C.5*H 
GO TO 2 

THIS RECYCLES THE INTEGRATION IF THE TRUNCATION ERROR IS EXCESSIVE. 

11 IF ( U2 - 0.031 ) 12, 13, 13 

12 DT = 2.0*H 
GO TO 14 

13 DT = SQRTFt SCRTF (0.5/U2) )*b 

14 IF (DT - A( ICO) <0.001 ) 16 ,16,15 

15 DT = A ( 100) *C.CC 1 

THIS SETS THE MAXIMUM STEP SIZE TO 1.0E-03 TIMES THE TOTAL TIME. 

16 T = T + 2 . 0* H 
C ( 1C4 ) = DT 
DO 17 I = 1 , N 

17 X ( I ) = XS ( I ) 

THIS UPDATES T AND X ( I ) . DT IS UPDATED BY STATEMENT 8,12,13 OR 15. 

RETURN 

END 

FUNCTION RELAY (R,DZONE,V) 

I F ( A BSF ( R ) - DZCNE)3,3,4 

3 RELAY = 0.0 
RETURN 

4 RELAY = S I GNF ( V , R ) 

RETURN 

END 

SUBROUTINE DERIV (T,X,XDCT,C) 

DIMENSION X ( 150) ,XDOT( 10 0),C( 150) 



noon ooo 



Fig. 40. Subroutine Deriv for Program Analog 



COMMENTS 

SUBROUTINE FCR CPEN - CLOSED CYCLE SYSTEM USING A TIMING C K T WHERE 
TIME TO SWITCH AND TIME TC TURN OFF IS PROP TC THE COMMAND SIGNAL 
AND THE RELAY OUTPUT IS A CONSTANT IN THE NONLINEAR REGION 
R = C ( 1 ) 

ERROR = R - X ( 1 ) 

IF(SWl) 599, 599, 600 



599 


IF { ABS F ( ERROR - C(2) ) ) 


603, 


600 , 


600 


600 


SIGN = 


ABSF ( ERROR ) /ERROR 










SW1 = 


1.0 










SW2 = 


O.C 










IF ( T 


- TS - TIME) 601 , 


602, 


602 




601 


DRIVE 


= C ( 6 ) * SW1 










GO TO 


604 








602 


DRIVE 


= -C { 6 )* SW1 










I F ( ABSF ( XDOT ( 1 ) ) - .05) 


603, 


603, 


604 


603 


SW1 = 


0.0 










S W2 = 


1.0 










TIME = 


X(0) 








604 


POWER 


= C ( 3 ) * ERROR * SW2 + 


DRIVE 


* SWl 



XDOT (1 ) = X ( 2 ) 

XDOT ( 2 ) = POWER * C(4) - X { 2 ) * C ( 5 ) 

TS = 0. 075 * ABS F ( R ) + 0 .05 
TO = 0.092 * ABSF(R) + 0.105 
X ( 3 ) = SW1 

X ( 4 ) = SW2 
X ( 5 ) = DRIVE 
X ( 6 ) = POWER 
X(lOl) = ERROR 
X( 102) = —XDOT ( 1 ) 

RETURN 



0 MM ENTS THIS IS BLUE DECK ONE 

' X ( 1 ) =SYSTEM OUTPUT , X(2)=SYSTEM RATE, C(l)= INPUT, 
- OF LINEAR ZONE, C(3)=SYSTEM GAIN, C(4)=PLANT GAIN, 
C ( 6 ) =MA GN I TUCE OF DRIVING VOLTAGE IN NONLIN REGION 
END 



C ( 2 ) = 1/2WIDTH 
C { 5 ) = 1 /TIME C 



Fig* 4l- Data input for Program Analog 



Kill tKMAN NALLE Y 
02 EQUATIONS 
04 RUNS 

01 TITLE CARDS 
OPEN-CLOSED LOOP 
ZERO COEFFS 
06 COEFF. CARCS 



IHES I S 



SERVC SYSTEM WITH SWITCHING = F ( T I ME , 



1 

2 

3 

. 4 

5 

6 

ZERO 
00 IC 
READ 
0.0 



0.5 
0.25 

40.00 

1 .0 

5.0 

100.0 

ICS 

CARDS 
TIME DATA 
1 . 5 



COMPUTE STEP SIZE 

-1 IS ORDER OF SMALLEST VARIABLE 
READ PRINT CATA 
07 VARIABLES PRINTED 

0 1 2 3 4 5 6 7 

02 TITLE CARDS 

TIME SYSTEM SYSTEM SWITCH SWITCH 
OUTPUT RATE ONE TWO 

20 INCREMENTS BETWEEN PRINTS 
READ GRAPH DATA 
02 GRAPHS 

1 0 102 101 

01 INCREMENTS BETWEEN POINTS 
HOLD COEFFS 
01 COEFF. CARDS 
1 1.0 
ZERO ICS 

00 IC CARDS 
HOLD TIME DATA 
HOLD STEP SIZE 
HOLD PRINT DATA 
HOLD GRAPH DATA 
HOLD COEFFS 

01 COEFF. CARDS 

1 1.5 

ZERO ICS 

00 IC CARCS 
HOLD TIME DATA 
HOLD STEP SIZE 
HOLD PRINT CATA 
HOLD GRAPH DATA 
HOLD COEFFS 

01 COEFF. CARCS 

1 2 .0 

ZERO ICS 
00 IC CARDS 
HOLD TIME CATA 
HOLD STEP SIZE 
HOLD PRINT DATA 
HOLD GRAPH DATA 



CRI VE 



PCWER 



- 72 - 



INPUT ) 



Fig. 42. Typical output for Program Analog 



OPEN-CLOSED LOOP SERVO SYSTEM WITH SWITCHING = F ( T I ME » INPUT) 
4 RUNS ARE CALLED FOR. 



RUN NUMBER 1 

THE NON-ZERO DATA COEFFICIENTS ARE AS FOLLOWS 



C( 


1 ) 


= 


. 50000 E+OD 


C( 


2) 


= 


.25000E+00 


C( 


3) 


S 


.4 0000E + 02 


c< 


4) 


— 


. 1 0000E+01 


C( 


5) 


= 


. 50000 E + 01 


C( 


6) 


= 


. I0000E+03 



THE NON-ZERO INITITAL CONDITIONS ARE AS FOLLOWS 
NONE 



THE TIMING DATA ARE AS FOLLOWS 

INITIAL TIME = .OOOOOE+DO 
FINAL TIME = .15000E+D1 

THE STEP SIZE IS COMPUTED, BASED ON 
SMALLEST VARIABLE OF ORDER l.OE-1 



PRINT SUMMARY 

5 INCREMENTS BETWEEN PRINTOUTS 
THE VARIABLES PRINTED ARE 



X( 


0) 


X( 


1 ) 


X( 


2) 


X( 


3) 


X( 


4 ) 


X( 


5) 


X( 


6) 



* X(O) REPRESENTS TIME 

GRAPH SUMMARY 

NO GRAPHS 
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KITTERMAN MALLEY THESIS 
RUN NUMBER 1 



TIME 



.OOOOOE+OO 
. 93000E- 03 
.128108-01 
.278 1 0E-01 
. 42 8 1 OE-0 1 
.578 1 OE-O 1 
.7281 OE-O I 
.8781 0E-01 
. 10281E + 00 
. 1 1781E+00 



. 1 3281 E+00 
. 1 4500E + 00 
. 14500E+00 
. 1 4500E+00 
. 14500E+00 
. 14502E + 00 
. 14579E+00 
. 1 57 1 5E + 00 
. 1 72 1 5E+00 
.■18715E + 00 



.202 1 5E + 00 
.2171 5E+00 
. 232 1 5E+00 
.2471 5E+00 
.262158+00 
.2771 5E + 00 
. 292 1 5E+00 
. 307 1 5E + 00 
.3221 5E+00 
.3371 5E + 0C 



.35215E-S-00 
.3671 5E+00 
. 382 1 5E + 00 
. 397 1 5E+ 00 
.4121 5E + 00 
.427 1 5E+00 
.4421 5E+00 
. 4571 5E+00 
.4721 5E+00 
.4871 5E+00 



. 502 1 5E + 00 
. 5 1 7 1 5E + 00 
.53215E+00 
.54715E+00 
. 562 1 5E + 00 
.5771 5E+00 
. 5921 5E + 00 
.60715E+00 
.6221 5E + 00 
.6371 5E+00 



SYSTEM 

OUTPUT 



.OOOOOE+OO 
-.43173E-04 
-. 80324E-02 
-. 36938E-01 
-. 85432E-01 
-. 1521 OE+OO 
-. 23563E+00 
-. 33479E+00 
-. 44847E+00 
-. 57562E+00 



-. 71 525E+00 
-. 83727E+00 
- . 8 3730E +00 
83730E+00 
- . 8373 1 E+00 
-.83755E+00 
-. 84538E+00 
- . 95485E + 00 
-.10790E+01 
-. 1 1805E+01 



12601E+01 
-. 13187E+01 
-. 13574E+01 
-. 1 3771 E + 0 1 
-. 1 3791 E+01 
-. 13648E+01 
-. 1 3353E+0 1 
1 2920E+01 
-. 12363E+01 
-. 1 1696E+01 



-. 10933E+01 
-. 10087E+01 
-.91707E+00 
-. 81 983E+00 
-.71817E+00 
-.61331 E+00 
-.50638E+00 
-. 39845E+00 
- . 2 9054 E+ 00 
-. 18358E+00 



-.784228-01 
.241498-01 
. 1 2343E+00 
. 2 1 880E+00 
. 30972E+00 
. 3957 1 E+00 
. 476 40 E + 00 
. 551 46E+00 
. 62065E+00 
.683808+00 



SYST EM 
RATE 



.OOOOOE+OO 
-.927 84E-0 1 
-. 12408E+01 
- . 259 63E +0 1 
-.38538E+01 
- . 502 05 E + 0 1 
-.6 1029E+01 
-. 7 1 0 70E + 01 
-.803 86E + 01 
-. 89029E+0 1 



- . 970 48E + 0 1 
-. 10313E+02 
-.103 1 4E+02 
-.1031 4E+02 
-.103 1 3E+02 
-.10311 E+02 
-. 10231 E+C2 
- . 904 69E +0 1 
-.751 48E+01 
-. 602 82 E+0 1 



-. 45969E+01 
-.32291 E+01 
-. 19322E+01 
-.712 32E + 00 
.425 63E+00 
.14777E+01 
.244098+01 
.331 34E + 01 
. 409 4 1 E + 0 1 
. 473 29E + 0 1 



. 533 05E + 01 
.583 84E + C 1 
. 6 30 85E + 0 1 
. 664 36E + 01 
.689 69E+0 1 
.707 22 E+0 1 
.717 35E+0 1 
.720 53E + 01 
.71725E+01 
.707 98E + 01 



. 693 25E + 01 
. 673 58E + 0 1 
.649 49 E + 0 1 
.621 52E+01 
.590 18E+01 
. 553 99E + 0 1 
.519 46 E+01 
.481 07E+01 
.441 28E + 01 
.400 54E + 0 1 



SWITCH 

ONE 



. OOOOOE+OO 
. 1 OOOOE+O 1 
. 1 OOOOE+O 1 
. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 OOOOE + O 1 
. 10000E+01 
. 10000E+01 



. 1 OOOOE+O 1 
. 1 0000E+01 
. 1 OOCO E+0 1 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 



.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOCOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 



.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 



.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
•.OOOOOE + OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOOOE+OO 



SWITCH 

TWO 



. 10000E + 01 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. 00000c - 90 
. OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 



.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. 1 0000 E + 0 1 
. 10000E+01 
. 1 OOCOE+O 1 
. 1 OOOOE+O 1 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 OOOOE+O 1 



. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
. 10000E + 01 
. 1 OOOOE+O 1 
. 1 0 000 E+ 0 1 
. 1 0000 E+0 1 
. 10000E+01 
. 1 OOOOE+O 1 



. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 0000 E + 0 1 
. 10000E+01 
. 1 0000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 OOOOE+O 1 



. 10000E+01 
. 10000E+01 
. 1 OOOOE+Ol 
. 1 OOOOE + Ol 
. 1 OOOOE+O 1 
. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 1 OOOOE+Ol 
. 10000E+01 



DR I VE 



. 1 OOOOE + O 3 
. 10000E+03 
. 1 OOOOE + O 3 
. 1 0000E+03 
. 10000E+03 
. 10000E + 03 
. 1 00005+03 
. 1 OOOOE + O 3 
, 1 0000E+03 
. 1 0000E + 03 



. 1 0000 E + 0 3 
. 1 0000E + 03 
. 10000E + 03 
. 1 0000 E+0 3 
. 1 OOOOE+03 
. 1 0000E+03 
. 1 OOOOE+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 



. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE + O 3 
. 1 OOOOE + O 3 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE + O 3 
. 1 OOOOE + O 3 
. 1 0000E+03 



. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE + C r 
. 1 0000E+03 
. 1 OOOOE+O 3 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE+O 3 
. 1 0000E+03 
. 1 0000E+03 



. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE + O 3 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE + O 3 
. 1 OOOOE+O 3 
. 1 0000E+03 
. 1 OOOOE+03 



i. 
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